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The Einstein - de Broglie particle-soliton concept is applied to simulate stationary states 
of an electron in a hydrogen atom. According to this concept, the electron is described 
by the localized regular solutions to some nonlinear equations. In the framework of Synge 
model for interacting scalar and electromagnetic fields a system of integral equations has 
been obtained, which describes the interaction between charged 3D soliton and Coulomb 
center. The asymptotic expressions for physical fields, describing soliton moving around 
the fixed Coulomb center, have been obtained with the help of integral equations. It 
is shown that the electron-soliton center travels along some stationary orbit around the 
Coulomb center. The electromagnetic radiation is absent as the Poynting vector has non- 
wave asymptote 0{r~^) after averaging over angles, i.e. the existence of spherical surface 
corresponding to null Poynting vector stream, has been proved. Vector lines for Poynting 
vector are constructed in asymptotical area. 
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1 Introduction 



From the history of quantum mechanics it is known that as early as 1927 in the framework 
of his "theory of double solution" Louis de Broglie made an attempt to represent the 
electron as a source of waves obeying the Schrodinger equation [1]. Later he modified 
his model showing that the electron should be described by regular solutions to some 
nonlinear equation coinciding with the Schrodinger one in the linear approximation. This 
scheme became famous as a causal nonlinear interpretation of quantum mechanics [2]. 
Developing this concept, de Broglie remarked that it had much in common with Einstein's 
ideas about unified field theory according to which particles were to be considered as clots 
of some material fields obeying the nonlinear field equations [3]. In recent years, these 
types of field configurations, known as soliton or particle-like solutions, came into active 
use to model extended elementary particles [4]. 

In this paper the Einstein-de Broglie soliton concept is employed to model stationary 
states of the electron in a hydrogen atom. 

2 Bohm problem about nonlinear resonance and its 
possible solution 

As a starting point, we will consider an interesting problem posed by D. Bohm. Long 
ago, in his book [5] D. Bohm discussed the possible connection between the wave-particle 
dualism in quantum mechanics and the hypothetical nonlinear origin of fundamental 
equations in a future theory of elementary particles. To illustrate the line of Bohm's 
argument we will consider a simple scalar model in the Minkowski space-time given by 
the Lagrangian density 

£ = di(i)* dj(i)rf^ - {mc/hf(j)*(j) + F{(j)*(j)). (2.1) 

Here i,j — 0, 1, 2, 3; 77'^ = diag(l, —1, —1, —1), the nonlinear function F{s) behaves at 
s ^ as s"^, n > 1, and is assumed as such that the corresponding field equations allow 
the existence of particle-like (soliton) solutions, i.e. regular configurations localized in 
space and endowed with finite energy. In particular, it can be shown that if one chooses 
F{s) — ks^/'^, k > 0, the model (2.1), known as the Synge one [6], admits the following 
stationary solutions: 

4>o — u{r) exp {—iujot) , r = |r|. (2.2) 

Here, the real radial function u{r) is regular everywhere and exponentially decreases as 
r — > 00, that provides finiteness of energy of the configuration 

Jd'xT'\<Po), (2.3) 

where T*-^ is the corresponding energy-momentum tensor. 

Moreover the model mentioned is intriguing due to the fact that nodeless solitons 
turn out to be stable by Lyapunov provided that their charge is fixed [7]. So there exist 
perturbed solitons slightly differing from the stationary sohtons (2.2): 

(/. = 0o + e(t,r). (2.4) 

Note that the perturbation ^ in (2.4) is small as compared with 0o only in the area of 
localization of the soliton, where 0o significantly differs from zero. None the less far from 
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the soliton center, where 0o is neghgibly small, one can put ~ ^, i.e. the tail of the 
soliton is completely defined by the perturbation 

D. Bohm put the following question: Does there exist any nonlinear model for which 
the spatial asymptote (as r — > oo) of a perturbed soliton-like solution represents oscilla- 
tions with characteristic frequency uj = E/W. In other words, for the model in question 
the principal Fourier- amplitude in the expansion of the field ~ ^ as r — > oo should 
correspond to the frequency uj connected with the soliton energy (2.3) by the Planck-de 
Broglie formula 

E = huj. (2.5) 

Note that for the model (2.1) at spatial infinity, where 0^0, the field equation 
reduces to the linear Klein- Gordon one 

[□ - {mc/hf](t) = 0, (2.6) 

and therefore the relation (2.5) holds only for solitons with unique energy E = mc^ defined 
by the mass m fixed in (2.1). Thus the universality of the relation (2.5) breaks down in 
the model (2.1), so forcing its modification. In the light of the above universality, the 
frequency uj in (2.5) being defined by the mass of the system, it seems natural that in the 
new, modified model one should use the gravitational field, spatial asymptote of which 
is also defined by the mass of the considered localized system. Thus, to solve the Bohm 
problem the possibility to invoke the gravitational field comes into reality [8] . 

So we will describe the new model with the Lagrangian density C — + Cg, where 

Cg = c^R/IGttG 

corresponds to the Einstein's theory of gravity, and Cyn is chosen as 

jCm = di<f)* d,<Pg'' - I{g^jWct> + F(0*0). (2.7) 

The crucial point of this scheme is to build up the invariant I{gij) depending on the metric 
tensor Qij of the Riemannian space-time and its derivatives. This invariant should possess 
such properties that in the vicinity of the soliton with mass m, the relation 

]im I {g,j) = {mc/hf (2.8) 

should hold. It is easy to see that on the basis of (2.8) one can asymptotically deduce the 
equation (2.6) from the Lagrangian (2.7). 

We argue that the invariant / can be built from the curvature tensor Rijki and its 
covariant derivatives Rijki-n' 

I = {I^Jll)c^h-^G-\ (2.9) 
where G is the Newton's gravitational constant and the invariants /i and I2 take the form 

h = RijkiR'''''/^S, h = -i?^,w;ni?'^'''7432. 

Estimating R'-^^^ at large distance r with the help of the Schwarzschild metric, one can 
find 

7i = GW/(cV«), 72 = GW/(cV«). 

So from (2.9) there immediately follows (2.8). Thus, within the modified model (2.7) 
for all massive particles the Planck-de Broglie relation (2.5) is automatically fulfilled. 
It means that in the framework of the scheme mentioned the principle of wave-particle 
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dualism is valid, according to which the relation (2.5) is realized as a condition of the 
nonlinear resonance. 

To verify the fact that solitons can really possess wave properties, the gedanken diffrac- 
tion experiment with individual electron-solitons similar to the numerical one of Biberman 
et al. [9], was realized. Solitons with some velocity were dropped into a rectilinear slit, 
cut in the impermeable screen, and the transverse momentum was calculated which they 
gained while passing the slit the width of which significantly exceeded the size of the 
soliton. As a result, the picture of distribution of the centers of scattered solitons was 
restored on the registration screen, by considering their initial distribution to be uniform 
over the transverse coordinate. It was clarified that though the center of each soliton fell 
into a definite place of the registration screen (depending on the point of crossing of the 
slit and the initial soliton profile), the statistical picture in many ways was similar to the 
well-known diffraction distribution in optics, i.e. Fresnel's picture at short distances from 
the slit and Fraunhofer's picture at large distances [10,11]. 

Fulfillment of the quantum mechanics correspondence principle for the Einstein-de 
Broglie's soliton model was discussed in the works [12-15]. In these papers it was shown 
that in the framework of the soliton model all quantum postulates were regained at the 
limit of point particles so that from the physical fields one can build the amplitude of 
probability and the average can be calculated as a scalar product in the Hilbert space by 
introducing the corresponding quantum operators. In this paper, we will show that in the 
framework of the Einstein-de Broglie soliton model a hydrogen atom can be simulated. 



3 Fundamental equations and structure of solutions 

z 



Let us consider the hydro- 
gen atom with the electron 
replaced by a localized ob- 
ject "soliton" that is mov- 
ing round the nucleus. So 
that the soliton-like solution 
does exist one has to con- 
sider a nonlinear model. 




As physical fields we choose the complex scalar field (p interacting with the electro- 
magnetic one = diAk — dkAi. The nucleus field is assumed to be the Coulomb one: 
A?^* = S'^Ze/r. The Lagrangian density is taken in the following form 

C = --^{F,k)' + \[dk- ie{Ak + AfM' - {mc/hfr^ + F{ct>*ct>), (3.1) 
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where e = e/{hc) is the couphng constant, is some nonhnear function, decreasing 

faster as ^ than |0p and is chosen so that the field equation at = allows the 
existence of stable stationary soliton-like solutions of the type (2.2), describing configu- 
rations with mass m and charge e. 

Note that for simplicity we do not write down the terms corresponding to the grav- 
itational field that will be taken into account implicitly with the help of the nonlinear 
resonance condition (2.5). 

Let us consider the nonrelativistic approximation assuming that 

= V exp {-imcH/n), (3.2) 

neglecting in the equations of motion higher derivatives of ip with respect to time and 
retaining only linear terms in A^. As a result, taking (3.2) into account we get the following 
system of equations: 

ih dti^ + {h'^/2m)Ax/j + (ZeVr)V' = -(^V2m)/(A, Ao, (3.3) 

aAo = {87^me/h^)\^l)\'^ = -Airp, (3.4) 

□A = 47r[2e^A\ilj\'^ - ie{ilj*Vi/j - i/jVi/j*)] = -(47r/c) j, (3.5) 

where 

f(A,Ao,ijj*ijj)iP = 2ie(AV)V' + 2(emc/;i)Ao'0 + ie'0divA + F'(V'*'0)V'- 

Moreover in the equations (3.3 - 3.5) it is supposed that the 4-potential Ai of the proper 
electromagnetic field of the soliton obeys to the Lorentz condition 

dtAo + c divA = 0, 

which is consistent with equations (3.3) - (3.5) owing to the conservation of electric charge. 

We will seek for the solutions to equations (3.3) - (3.5) describing the stationary state of 
an atom when the electron-soliton center is assumed to be moving along a circular orbit of 
the radius oq with some angular velocity Q. In this problem there arise two characteristic 
lengths: the size of the soliton / = h/{mc) and the Bohr's radius a = / [mZe^]. It is 
obvious that oq ~ a ^ 

Let us first consider the area near the soliton center where r — oq ~ Suppose the 
soliton center trajectory to be r = C,{t). Putting into (3.3) the configuration 

■0 — u{y — C(^)) exp (iS/h), 

neglecting the contribution of the proper electromagnetic field and separating the real 
and imaginary parts, we get 

AS + 2(VS-mC)-Vu/u^0. (3.7) 

Assuming 5* to be a slowly varying function of a point in the vicinity of the soliton center, 
from (3.7) we deduce 

S ^mC ■ {r - C) + Cot + x{t), Cq = const. (3.8) 
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Taking into account the classical equations of motion of a charged particle in the Coulomb 
field 

and using the expansion 

1^.1 C-(r-C) 

from (3.6) and (3.8) we derive 

5a = yC + ^ = /:(t), 

where L{^) is the Lagrangian of a particle in the Coulomb field. Thus, the function x is 
the classical action on the trajectory: 



X(t) = I C{t)dt, (3.9) 





and the function u is the soliton-like solution to the quasi-stationary problem 

n'^{f + Au/u) = 2mCo. (3.10) 
In this case according to (3.4) and (3.5) 

p = -{2me/h'^)u^, j = -2ecM^(eA + mC/h), 

which makes it possible, using the common solutions to equations (3.4), (3.5) and (3.10), 
to calculate the potentials of the electromagnetic field in the vicinity of the soliton 
center: 

Ao = Ao{r - C(t)), cA = C(t)Ao(r - C(t)), 

where the terms C /c^ are neglected. 

Let us now study the asymptotic behavior of soliton at large distance, i.e. we will 
study it's "tail". We will use successive approximation method. Thus, we need to rewrite 
the differential equations (3.3), (3.4) and (3.5) in integral form. Note that we are not 
solving Cauchy evolution problem, choosing definite initial condition. Were it a question 
of Cauchy problem, we would bound to use retarded Green's function. In this case it has 
been assumed that the object under consideration (atom) already exists infinitely long. 
Thus we consider the problem to study corresponding steady state, which ehminates the 
possibility of using retarded electromagnetic Green's function. The problem mentioned, 
according to our view, can be satisfied by half-sum of retarded and advanced solutions. 
Above mentioned selection can be justified by the assumption that during the evolution 
process the radiation of the independent i.e. half- difference of retarded and advanced 
electromagnetic fields took place. 

To find out the field ip far from the soliton center, we rewrite equation (3.3) in the 
integral form 

ip{t,r) = Cnipnir) cxp {-lUJ J) + 

+ ^J^^J(^^'J exp[-iuj{t - t')]G{r,r';uj + iO)fip{t' ,r'), (3.11) 
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where V'n(r) is the eigenfunction of the Hamiltonian of a hydrogen atom for a stationary 
state of number n with energy En — hcUn, Cn — const and G(r, r'; uj) is the Hamiltonian's 
resolvent having the form [16] 



G(r,r ;w) - 



Wii,^i/2{-ikr+) Mjj.,i/2(-iA;r_) 



(3.12) 



AttR 

Here,the following notation is used: 

k = {2mou/ny/^, lmk>0, i/ = (ka)'^, 
r± = r + r ±|r — r| 

and the Whittaker functions H^iiy,i/2, -^ii/,1/2 and their derivatives W^ji/,1/2, Mii,^i/2 are in- 
troduced. To find the field at large distances from the electron-soliton center, i.e. at 
|r — ao| ^ l-i it is sufficient to put in (3.11) 

/V^(t,r) =51 exp(-ia;„t)5(r-C(t)), g = const, (3.13) 

where the relation (2.5) is taken into account. As a result, we get 

ii{t,Y) = C„V^„(r) exp (-io'nt) + 

+ — j duj j dt' exp [-iujt + it'{uj - uj^)] G{r, r';u; + iO). (3.14) 

It is easy to verify that the field (3.14) decreases exponentially at large distances. With 
the help of (3.14) and equations (3.4) and (3.5), one can evaluate the electromagnetic 
field outside the soliton. In (3.14) C„ is unknown constant and ■0n(r) is the wave function 
of electron in stationary state. In each step of iteration one obtains (3.14) where the 
stationary "tail" of soliton is marked out and its center moves along some effective orbit. 
The orbital parameter and constant C„ may be defined in arbitrary approximation of 
minimization | [■^(j^) —ip(^k^i) \ \ [17] . Still now the constants C„ and g are not found explicitly 
and we hope to obtain them in our forthcoming papers. 

Let us now solve the equations (3.4) and (3.5). Considering that the nonlinear source 
is rather weak one, we will replace the right hand sides of the equations (3.4) and (3.5) 
by 5- functions. Let us also notice that 

E_ = (E_ + E+)/2 + (E_ - E+)/2, 

where E_ = W^^, E_|_ = E"'^''. It is well known that the half-difference of retarded and 
advanced fields radiates. So it will be sufficient to consider, as was discussed earlier, the 
half-sum of retarded and advanced solutions, describing the stationary state. It means, 
we will seek the strength of the electromagnetic field as the half-sum of those for retarded 
and advanced fields. It means that for large times \uJn\t ^ 1 the 4-potential A/- will 
contain only stationary part = {A'^*- + A'ff^)/2. 

Let us find the expression for E_. Radius of the soliton I is rather small in comparison 
to Bohr radius a, i.e. a ^ I. So the source can be considered as proper one. Let the 
point-like charge e move along the given trajectory r = (^(t) with velocity v(t) = C(t)- 
Then, to describe the electromagnetic field, generated by the charge, one can write charge 
density and current density as 

p{t,r)^e5[r-Cit)], 3{t,r) ^ ey{t) 5[r - C{t)]. (3.15) 
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Then the equations (3.4) and (3.5) take form: 

□^0 = 
□A = 



47re5[r-C(t)], 



that lead to the weU-known Lienard- Wiechert potentials. As we know, the retarded time 
writes 

t-^t- R{t-,r)/c, 

which leads to 



1 



dt 1 — (n_ ■ v)/c' 



(9R(t_,r) 

ml 



= -(n_-v), Vi- = 



n_ 



c — n_ • V 



where n = Hit, r)/R{t, r). Later, using the following expressions for Ao_(i,r) and 
A_(t,r) [18,19] 



ec 



[cR - (v • R)]t. 
one finds the expression for retarded field 



-, A_{t,r) 



ev 



cR-{v- R) 



(c^ — f ^)(cn_ — v) ^ [n_, [(cn_ — v), v] 



and 



i?2[c- (n_ • v)]3 i?[c - (n_ • v)]3 

B_ = [n_,E4_. 



(3.16) 



(3.17) 



Analogically writing the advanced time as 

t+ = t + R{t+,r)/c, 

for the advanced field we find 

(c^ - v2)(cn_ + v) [n_, [(cn„ + v), v]^ 



E+ = e 



i?2[c+ (n_ ■ v)]= 



+ 



R[c+ (n_ ■ v)]' 



and 



B_ 



■[n+,EH 



(3.18) 



(3.19) 



To calculate the power, lost by the charge due to radiation, one has to compose Poynting 
vector S and retain the terms of the order 1/i?^ as the integration will take place along 
infinitely distant surface. As we know, Poynting vector is expressed by the relation 



S ^ -[E,B1. 



In this case the field strengths are 



E = i(E+ + E_), B = ^([n_, E_] - [n+, E+]), 



(3.20) 



(3.21) 



where E_ = E''^*, E+ = E"'^", n± = R±/-R±, R± = r - C(^±) and t± are the roots 
of the equations t± — t ± R±/c. 
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Using some manipulations from vector analysis we rewrite Poynting vector as 



S = T^{n-Ei - n+E^ + (n_ - n+)(E_ • E+) - 

-E_((E_ + E+) • n_) + E+((E_ + E+) • n+)}. (3.22) 

Let us now rewrite all these in spherical system of coordinates. In this case 

R(r, t-t) = {r — flo sin^C0SQ;±, — Oo cos^cosq;±, Oo sina-t}, 
= {ao ^ sinO sma±, aofl cos9 sina-i-, aoflcosa±}, 
v(t-i-) = {— ao sin6' cosq;±, — ao cos^cosq;±, Oo sincci}, 



R{r,t±) — ^Jr'^ + Oq — 2aor sin^cosQ;±, 



where a± — a — D,t±. Retaining the terms of the order oo/r for n_ and n+ one can 
obtain 

n-t = {1, — (ao/r) cos6' cos(q; — flt±), (ao/r) sin(Q; — flt±)}. 

Using the first nonvanishing approximation of the order v/c, from (3.16) and (3.18) we 
will get 



E_ fa 
E+ f=:i e 

Taking into account that 



n_ n_(ii • v) v 



n+(n+- v) -V 



i?2 cm 



(3.23) 
(3.24) 



1/R ^ 1/r - (ao/r2)sin^cosQ;_, 1/R^ ^ 1/r^ 
(n_ • v) -aofi^ sin6' cosq;_ + {alQ'^/r)[l - sin^6'cos^Q;_], 



where q;_ = a — for E_ one gets 

1 
r 

Analogically one finds 



E_ f« e{ — [1 + {alH'^/c^){l - sin^6'cos^Q!_)], (ao^^/c^r) cos^cosq;_, -{aoft^/c^r) sina_}. 



E+ fti e{ — [1 + {alQ,'^/c^){l - sin^6'cos^Q!+)], (ao^^/c^r) cos6'cosq;+, -{agO,'^ /c^r) sina+j. 

Putting the above expressions for n±, E± into (3.22) one can find the expressions for 
Poynting vector. In doing so we will take into account that the normals n_ and n_|_ 
coincide as r — ^ oo with n = r/r. Retaining the terms {ao/r)^ and also taking into 
account that agQ = v <^ c, for the circular motion in the spherical coordinates r, 6 and a 
we have the following components of the Poynting vector S: 

g2 ^2 ^4 

Sr = V^sin2^sin2(Q;-Qt)sin(2fir/c), 

Se — — ^ cos^sin(Q; — Qi) sin(Jlr/c), 

Sa = , ° , cos(a -Ot) sin (Or/c). (3.25) 
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It is obvious that for = ckn/Q. with k — 0, 1, 2, • • • all the components of Poynting 
vector turn to zero i.e. S = 0. 

All the calculations done above can be summed up as follows. From (3.21) it follows 
that the projection of the Poynting vector S in the direction of the vector N = (n_|_+n_)/2, 
coinciding as r ^ oo with n — r/r, takes the form 

•^iv = - + n+ • n_)V2. (3.26) 

Since n± = n + 0{r~^), after averaging expression (3.26) over the sphere, we find 

< >^ ^(^< E'_> - < El>~^ ^ 0{r-'). (3.27) 

Thus according to (3.27) the electromagnetic radiation from the system is absent. In 
particular, for the circular motion in the spherical coordinates r, 6 and a we have the 
following structure of the Poynting vector S: 

Sr = ^sm'^9sm2{a-nt)sm{2nr/c), (3.28) 
So = sin(Qr/c)0(r-^), -S^ = sin(Qr/c) ©(r^^), 

where k = e'^alfl'^/{16Tic^). From (3.28) as well as from (3.25) it is obvious that there 
exist spherical surfaces where either 5"^ = or S = 0, thus once again confirming the fact 
that in the stationary states described, radiation is absent [20]. 

Let us describe the vector lines for Poynting vector. In spherical system of coordinates 
we have 

dr dO . ^ da 

= r- = rsm^— . (3.29) 

Last two fractions form integrable combination. Putting Sq and 5'a, from this equality we 
obtain 

tana da, (3.30) 



sin6' cos6' 
which leads to the first integral 

|tan^cos(Q; — ni)| = Pi, Pi = const. (3.31) 

Now we consider first two fractions. Putting cos(q; — Q,t) — Pi/tan^ we find 

dr 



Q? r^cos(r2r/c) 



= Pi sin^ de. (3.32) 



As we consider the region where r ^ cxd it is possible to factor out the term 1/r^ from 
the integrand. Then we obtain appoximately 

that leads to 

-— r ln|tan(7r/4 + VLrl2c)\ = -Pi cos^ + P2, P2 = const. (3.34) 
Thus, we built vector lines for Poynting vector in asymptotic region. 
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4 Conclusion 



In the considered soliton model of a hydrogen atom the stabihty condition of spatial 
stationary motions of electrons in the field of the Coulomb center is fulfilled. The existence 
of this kind of motion had also been anticipated by Boguslavsky [21] and Chetaev [22]. 
In particular, due to the fulfillment of the nonlinear resonance condition (2.5) the energy 
spectrum of these stationary states coincides with that of a hydrogen atom. This fact 
indicates the role of nonlinearity in the formation of extended micro-objects, whose laws 
of evolution agree with quantum mechanics. 
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